Transport quantum logic gates for trapped ions 

D. Leibfried, E. Knill, C. Ospelkaus and D. J. Wineland 

National Institute of Standards and Technology, 
325 Broadway, Boulder, CO 80305, USA 
(Dated: February 1, 2008) 

Abstract 

Many efforts are currently underway to build a device capable of large scale quantum information 
processing (QIP). Whereas QIP has been demonstrated for a few qubits in several systems, many 
technical difficulties must be overcome in order to construct a large-scale device. In one proposal 
for large-scale QIP, trapped ions are manipulated by precisely controlled light pulses and moved 
through and stored in multizone trap arrays. The technical overhead necessary to precisely control 
both the ion geometrical configurations and the laser interactions is demanding. Here we propose 
methods that significantly reduce the overhead on laser beam control for performing single and 
multiple qubit operations on trapped ions. We show how a universal set of operations can be 
implemented by controlled transport of ions through stationary laser beams. At the same time, 
each laser beam can be used to perform many operations in parallel, potentially reducing the 
total laser power necessary to carry out QIP tasks. The overall setup necessary for implementing 
transport gates is simpler than for gates executed on stationary ions. We also suggest a transport- 
based two-qubit gate scheme utilizing microfabricated permanent magnets that can be executed 
without laser light. 
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I. INTRODUCTION 



A number of physical implementations have been proposed for quantum information 
processing (QIP) This paper is based on a proposal where trapped atomic ion qubits 



are to be held in a large trap array M\- Its implementation requires transporting ions 
between separated zones, precise control of local potentials and, at the same time, precise 
control of laser beam pointing, intensity and pulse shape. These requirements create an 
imposing overhead of classical control for large trap arrays with multiple interaction zones. 
Ion transport is accomplished by electronically changing the potentials of individual control 
electrodes in the trap array [4] and might be realized with on-board CMOS electronics j^, 
a technology with a long and very successful track record for scaling. The situation is very 
different for the optics necessary for laser beam control: Microfabricated beam steering 
optics and electro-optical devices are typically still "one-of-a-kind" designs with only small 
numbers produced and scalability in the context of QIP still to be demonstrated. This 
problem is compounded by the wavelengths that are of interest in QIP with trapped ions, 
which are typically in the near UV between 214 nm and 400 nm. In addition, a mature 
optical fiber technology does not yet exist for this wavelength range. It is anticipated that 
high laser power will be required in QIP with trapped ions loj], so beam-splitters or lossy 
elements should be used sparingly. At the same time, fault tolerant architectures require 
implementing parallel operations 7|]. Miniaturization of the currently used approach with 
switched beams, as discussed in js] (for example), is based on the use of a large number of 
beamsplitters and control elements to achieve both parallel and individual control of many 
different gate operations. Even if the elements used have little loss (which is currently hard 
to achieve in the UV), parallel operations would magnify the already demanding power 
requirements. 

The purpose of this paper is to show that precise control of the time- dependent exter- 
nal potentials used to transport ions within a trap array can replace the requirement for 
precise temporal control of laser beam intensity in order to implement universal quantum 
computation. In an architecture based on transport, laser beams can be switched on and 
off collectively with relaxed requirements on timing and on/off ratios. Such a scenario may 
also allow for efficient use of one and the same laser beam in many parallel operations, thus 
achieving parallelism without the need for higher laser power. Under such circumstances 
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it is even conceivable to further enhance the available power in laser beams with optical 
cavities of modest finesse. 

The paper is organized as follows: Section [III outlines the basic architecture and QIP 
primitives necessary for universal quantum computation with the proposed scheme. We 
concentrate on qubits that are comprised of the hyperfine states of ions, which are manipu- 
lated by stimulated two-photon Raman transitions [2|], but it is possible to adapt the basic 
architecture for ion qubits of a different type. Section IIIII and IIVI briefly summarize the 
necessary Raman laser interactions and the spatial dependence of the laser beam modes 
used in subsequent sections. Sections |V] and IVH outline the details of one-qubit rotations 
and two-qubit phase gates implemented by ion transport through laser beams. Section IVIII 
discusses how sympathetic cooling can be incorporated and section IVIIII introduces some 
possible extensions of the scheme, including a two-qubit gate based on transporting ions 
over a periodic array of microfabricated permanent magnets, without the need for laser 
beams. Finally we summarize and offer some conclusions in section IIXI 

II. BASIC ARCHITECTURE 

The goal of the architecture discussed here is to minimize the requirements on laser-beam 
steering, pulse shaping and switching as much as possible by utilizing temporal control of 
potentials applied to the ions in a multizone trap array. Temporal control is already needed 
for efficient transport, separation and recombination of ions, so, with refinements, we can 
also employ it for qubit gate operations. 

Our proposed architecture is based on a multi-zone geometry [lOl. To be 



specific we consider planar surface electrode trap arrays [11] in the following, but the basic 
ideas should also work in other types of trap arrays. In this architecture, logic operations 
are implemented by two basic primitive steps: 

(i) The ions carrying the quantum information are arranged into a particular spatial 
configuration in the trap array while the laser beams are switched off [Fig. [1] (a)]. 

(ii) All laser-beam assisted operations scheduled for the configuration are implemented 
after (i) is carried out. The laser beam(s) are collectively switched on, then single qubit 
ions or pairs of qubit ions are transported through the laser beams to implement one-qubit 
rotations, two-qubit gates and measurements [Fig. [T] (b)]. Finally the laser beams are 
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switched off collectively. 

Steps (i) and (ii) are repeated until the computation is finished. In more detail, the 
control for ion motion in (i) can be accomplished with a few sequential elementary sub- 
steps. For example, these sub-steps could be translations of the potential wells containing 
an ion or ions in the array and splitting and recombing potential wells to reconfigure ions into 
different groupings 4, ll2|. These basic operations are indicated schematically in Fig. [T](a). 
In addition to these classical means of transport, quantum information can be transported in 
the array without physically moving the information carriers by teleportation 12, 1^, 14, 15|. 



Teleportation could be supported by a backbone of entangled qubits distributed over the 
whole array before and/or in parallel with the computation. Such an entanglement backbone 
could also be part of an efficient error correction scheme [l^. 

After the preconfiguration of qubits in the array, step (ii) is implemented. This step 
can be broken down into three basic laser assisted sub-operations that we call transport 
gates: single-qubit rotations, two-qubit gates and measurement. Single-qubit rotations can 
be initiated by first turning on specific laser beams globally over the entire array. Then 
the qubits scheduled for one qubit-rotations are transported through the laser beams in 
a controlled fashion (see Fig. [1] (b)). Next the beams for two-qubit gates are turned on 
and the pairs of ions scheduled for two-qubit gates are transported through the beams. 
Finally all qubits scheduled for measurement are read out by either turning on a detection 
beam at their current location or transporting them through a globally switched detection 
beam. Depending on the exact nature of the detection scheme, all measurements can be 
done in parallel if position-resolving detectors are used. Alternatively, detection could be 
accomplished serially with scheduled transports if no (or only limited) position resolution 
is available. The physical implementation of one-qubit rotations and two-qubit gates is 
discussed in sections |V] and |VT] respectively. 

Previously, temporal control of ions' internal states, such as qubit rotations, had to be 
achieved by individually calibrated, precise switching of laser beams. With transport gates 
laser beams can be switched on and off globally over the entire trap array while precise 
individual control is now transferred to the ion motion. This also facilitates the use of one 
set of laser beams for parallel operations on ions distributed over the trap array, reduces 
the complexity of optics, and might lead to lower requirements on the total laser power 
necessary to run processors with a certain number of qubits. 
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By repetition of (i) and (ii) we can realize one-qubit rotations, two-qubit gates and 
measurement between arbitrary qubits in an arbitrary order, which is sufficient for universal 
quantum computation jl^]. Individual operations are then controlled by the motion of ions 
alone, while the switching of lasers can be implemented with reduced timing precision. This 
could be of significant practical importance, if active feedback on the lasers to counteract 
intensity or beam position fluctuations is desired. In the traditional scheme where operations 
depend on the temporal characteristics of the laser beams, such feedback would have to act 
on timescales much shorter than that of the laser pulses. In the scheme discussed here, 
the transport of the ions can be delayed by a suitable amount of time for the feedback 
to settle. Such a procedure could also alleviate the detrimental effects of other switching 
imperfections such as phase chirps in acousto-optic modulators. As a further example of 
the potential simplification, consider the problem of implementing in parallel a specific 
rotation on several spatially separated qubits with the same laser beam. If the rotation is 
implemented by applying a pulse to ions already in place, we require the laser intensity to be 
the same at the site of each ion, a difficult task, given the general divergence/convergence 
of the laser beams. This problem can be alleviated by controlling the transport of each 
qubit through the beam separately. Furthermore, since the ions are transported completely 
through the laser beams, the gate interaction does not change if beams have small pointing 
instabilities in the plane of the trap array that change on a timescale long compared to 
the gate interactions (typically fulffiled for beam-steering timescales in the laboratory). It 
is therefore sufficient to stabilize the beam pointing in the direction perpendicular to the 
motion of the ions. 



III. TWO-PHOTON STIMULATED RAMAN INTERACTIONS 

In this section we briefly review the basic interactions that play a role in the transport 
gates discussed later. We assume that we have perfect control of the ion motion in the array; 
in particular, we assume that we can confine one or two qubit ions in one potential well with 
precisely defined motional frequencies, and we also assume we can translate the well along a 
predefined trajectory without changing the motional frequencies. We assume that the ions 
start in and remain in the motional ground state in the accelerating phases at the beginning 
and end of the transport and that the well is translated with a constant velocity v relative 
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to the laboratory frame while the ions move through the laser beams. This situation, as 
viewed from the ions' frame of reference, is equivalent to that of ions in the ground state in 
a stationary potential well under the influence of temporally controlled laser beams. We can 
therefore describe one-qubit gates and two-qubit gates with the interaction Hamiltonians 
that are also appropriate for ions at rest in the laboratory frame. The difference is that, 
for each operation, the coordinates of the ions are transformed by r — riab + v t with riab 
a point at rest in the laboratory frame that is taken to coincide with r at time t = 0. We 
assume the qubit state space is spanned by two hyperfine ground states of each ion, formally 
equivalent to the two states of a spin-1/2 magnetic moment in a magnetic field. We therefore 
label the qubit states |t) and |1) and denote their energy difference by hujQ . We consider 



stimulated Raman transitions [2, ll8| , implemented with two laser beams with wavevectors 
ki and k2 and frequencies ui and 002- After adiabatic elimination of the off- resonant levels, 
the effective interactions of one ion with two Raman laser beams can be described in the 
rotating wave approximation by the Hamiltonian j^, 19, 20 1 



Hj = hnoe-'^^^''>-^''^'+^^e'^^-''\ t)(i I +h.c. 

+ ;i^e-*(^"*+'^)e*^''""| t)(T I +h.c. 

+ hn^e~'^^'''~^^^e'^^-''\ I +h.c. 

+ ^AJi)(i|+^At|T)(TI, (1) 

where h.c. is the hermitian conjugate of the previous term, h is Planck's constant divided 
by 271, 60 = uji — UJ2 is the frequency difference, Ak = ki — k2 the wavevector difference and 
= 01 — 02 is the phase difference of the two laser fields at r = 0. The above expression 
breaks the total interaction down into the three parts: (i) interactions that can change the 
spin and possibly the motional state at the same time (sideband transitions) are associated 
with the Rabi frequency Qq, (ii) interactions that can change only the motional state are 
associated with the Rabi frequencies fig with se{], |}, and (iii) pure level shifts A^ due to 
the AC-Stark effect induced independently by each of the two laser beams. Because the 
energy of the ion's internal and/or external degrees of freedom changes in transitions of 
type (i) and (ii), these interactions can be only resonantly driven by one photon from each 
of the two electric fields (which also need to have an energy difference appropriate for energy 
conservation in the process) while the coupling in (iii) is mediated by two photons out of 
the same electric field (we do not consider the case (5o = 0). The Rabi frequencies are given 



by sums over dipole matrix elements: 



^ ^ (s|d-E2|0(^|d-Ei|s) , (s|d-Ei|/)(/|d-E2|s) 
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, _ 1 Ks|d-Ei|/)P |(.|d-E2|0P 
^ 4/1^ Y Aa A^2 ' ^ ^ 

where d is the dipole operator (assumed to be real by convention), Ei and E2 are the two 
(real) laser field amplitude vectors, and Aij = Uj — {Ei — Es)/Ti is the detuning of laser 
field i (je{l,2}) with respect to the near-resonant level |/). Typically | [) (| t)) is in the 
ground-state S manifold with energy {E^) while |/) is one of the levels in the first excited 
P manifold of the ion with energy Ei. In the regime of power and detunings of interest for 
this work [6|, the hierarchy ^ cuq ^ (fioj^s, A^) applies. In that case, the choices of 5q 
and Ak determine which terms dominate the evolution described by Eq. ([T]). For 5q ~ uoq 
the terms connected to process (i) are near- resonant, leading to spin-fiip transitions that 
can be highly independent of the motion if, in addition, the laser beams are co-propagating, 
leading to |Ak| ~ (see below). For 5q ^ and |Ak| ~ |k|, with ujm one of the motional 
frequencies of the ion(s), process (ii) dominates, leading to state-dependent driving of the 
motion. 

In both situations we would like to suppress (or at least precisely control) the AC-Stark 
shifts (iii) that are inevitable in the presence of the electric fields. In practice we can often 
balance the AC-Stark shifts A^ caused by process (iii) for | |) and | [) without diminishing 
the Rabi frequencies VLq and Vts appreciably by a judicious choice of the intensity and/or 
polarizations of Ei and E2. For example in ^Be"*" two polarizations close to linear and 
orthogonal to each other can be used to yield Aj — A| and Vti ~ — 2fi| for the states 
_[) = |F = l^mp = —1) and | |) = |F = 2, nip = —2) in the 5*1/2 electronic ground state 
20l |. On the other hand, we can use only one of the Raman beams (e.g. Ei 7^ 0, E2=0 with 
El circularly polarized) that induce state-dependent AC-Stark shifts to implement rotations 
around the z-axis of the Bloch sphere of the form (up to an irrelevant global phase) 

Z(0)|T) = e^<^|T) (3) 

m\i) = e-"^\i). 
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IV. PARAXIAL DESCRIPTION OF THE LASER MODES 



We assume that all laser beams can be described as having a Gaussian TEMqo transverse 
beam profile. A Gaussian beam with wavevector k and angular frequency uj propagating 
along the ^;-axis is the lowest order solution of the paraxial wave equation and can be 
described by [2l| 
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+ C.C. (4) 



(5) 



For simplicity we consider only the situation where \z\ -C \zr\, although it is possible to 
generalize the formalism to the curved wavefronts in the regions with \z\ > \zr\. In our 
region of interest we can write w{z) ~ wo, aj:ctaii{z / Zr) ~ 0, and {kr"^) / {2R{z)) ~ 0. Under 
these approximations the beam simplifies to 



E{r, ^) ~ ~^ exp[— r^/wg] exp {—i{kz — ujt + 0)) + c.c. 



(6) 



This expression describes the salient points of the spatial dependence of the beams in the 
following discussions. 



V. ONE-QUBIT ROTATIONS 

Up to the present, laser- induced one-qubit rotations of trapped ion qubits have been 
implemented by turning laser beams on and off for a duration appropriate to achieve a certain 
rotation angle on the Bloch sphere. By also controlling the relative phase of the laser beams, 
arbitrary rotations on the Bloch sphere can be implemented. Here we propose to move ions 
precisely through a laser beam to achieve the same level of control. The duration of the 
interaction is now controlled by the speed = |v| of the movement. We assume that the ion 
under study traverses two co-propagating laser beams at an angle 7 relative to the /c-vectors 
[see Fig. |2] (a)]. Since the wavevectors of the two laser fields are not exactly equal, their 
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relative phase changes by A0p ^ sq cuq/c = 27r(so/Ao) over locations separated by distance sq 
along the direction of the beams, where Aq is the wavelength corresponding to the hyperfine 
transition frequency of the ion in question. Typically ujq/(2'k) is one to several gigahertz so 
that in ^Be+, luq ~ 2n 1.25 GHz and Aq — 0.24 m. Therefore Ak ■ r = (ki — k2) ■ r can be 
taken to be constant for the typical variations of r of a few nanometers, and the effective 
Hamiltonian Eq. ([1]) is highly independent of the ion motion. If the laser frequencies differ 
by So = uJo, we can drop all rapidly rotating terms in Eq. ([1]) in a second rotating- wave 
approximation and are left with 

Hn^ = hnoe-"^\ DU l+h.c. (7) 

leading to Rabi rotations R{9, 0) on the Bloch sphere that can be described by 

Rie^m) = cos(^^/2)|T) -ze^^sin(^^/2)|i) (8) 
R{0,m) = -ze-^'^sin(e/2)|T) + cos(e/2)|i). 

For simplicity we assume that both beams have the same transverse mode function as 
given by Eq. ([6]). The trajectory is such that the ion is located in the center of the beams 
at t = 0, where it experiences the maximum coupling strength Qm, which can be calculated 
by inserting Eq of Eq. ([6]) into Qq of Eq. ([2]). The time dependent Raman- coupling strength 
is then 

n{t) = exp[-2(f sm{-f)t/wof]. (9) 

Integration of Eq. (Q gives the rotation angle 6 on the Bloch sphere resulting from a single 
traverse of the beams starting at time — T and ending at T: 

9 = r n{t) dt = ^]^^^^,/| eTi{V2vsm{^)T/wo), (10) 
J-T f sm(7) V 2 

where erf(z) = (2/y^) Jq e~^^dt is the error function. In practice, it is permissible to neglect 
effects due to the finite distance of the transport start and end points from the center of the 
modes (T oo), in which case erf (v^f sin(7)T/wo) — > 1 in Eq. ffTOj) . The relative error 
A9/6 introduced by this approximation is A9/6 = 1 — erf (-\/2-D sin(7)/wo)! where D is the 
initial distance of the ion from the center of the beams [see Fig. [2] (c)]. As an example, if 
7 = TT (beams perpendicular to the trajectory), the infidelity 1—F ~ {A6/6)^ is smaller than 
10~^ ii D > 2.6 wq. Based on this path length of 2.6 wq we can compare the approximate 
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duration Tt to execute a transport gate with angle 9 to the duration Ts necessary for the 
same rotation if the ion resides in the center of beams that are switched on and off suddenly, 



yielding Tt/T, = 2.6y'2/7r ~ 2.07. 

The exact value of 9 can be controlled by the speed with which the ion is transported 
through the beams. For example, to execute a vr-pulse given a waist size of = 20 /xm and 
f2m = 250 kHz, we need a velocity f ~ 25 m/s. Similar transport speeds (average speed 25 
m/s and a peak velocity of about 45 m/s), while keeping a single ion in the ground state of 
the transported potential well, have been demonstrated previously {4]. Moreover, the ground 
state does not need to be preserved for single qubit rotations (using co-propagating beams), 
so faster, nonadiabatic transport is possible if the ion can be sympathetically recooled after 
a rotation {3, 3|. 

The angle cf) on the Bloch sphere is controlled by the relative phase of the laser fields at 
riab- This phase changes by A0p ~ 27r(so/Ao) between two interaction zones utilizing the 
same beam that are separated by a beam path of length sq. Independent of whether or not 
temporal laser pulses or transport gates are used, A0p cannot be neglected over an extended 
trap array. It is necessary to do careful bookkeeping of phases over an extended trap array to 
properly set and control operations. If the same pair of beams is used for one-qubit rotations 
in several different places, we do not have the freedom of setting independently for all these 
rotations. There are several ways to deal with this problem. For example, a straightforward, 
but logistically costly, solution would be to transport a given ion to a location where (f) 
happens to be correct for the rotation scheduled on that ion. A more flexible strategy would 
be to have a discrete set of different relative phases that could be implemented by several 
beam pairs that are turned on in parallel to address separate locations or by sequentially 
switching the relative phase of the same beam (universal computation can be achieved with 
a discrete set of one-qubit rotations |22|). The time overhead caused by such sequential 
rotations may not slow down the computation considerably, because one-qubit rotations are 
typically much faster than two-qubit gates. 

We could also momentarily turn on only one of the Raman beams. Transporting the 
ion through that beam with a certain velocity induces a time-dependent AC Stark shift 
that can be used to achieve rotations Z{(j)) (see the discussion preceding Eq. ([3])). A given 
rotation R{92,4>) with (p fixed can be turned into an arbitrary one-qubit rotation U by 
sandwiching it between two Z-operations, U = e'^^'^Z{(j)^)R{92,(f))Z{(f)i) [Q]. In the context 
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of a multi-operation algorithm, this can be further simplified: The final operation Z{(j)3) does 
not affect measurement outcomes and can therefore be neglected in a one-qubit operation 
that precedes a measurement. This also holds if two-qubit p ha se gates are applied before 
the measurement since all Z{(j)) commute with phase gates [23]. For all other operations 
the final step of the previous rotation Zpi.ov(03) can be combined into the next one-qubit 
manipulation U' = where Z'{(p[) = Z(0;)Zprcv(03) = Z{^[ + ^3). Thus 

transporting the qubit ion through one laser beam first and then through two beams with 
detuning uq at the appropriate velocity is sufficient for implementing universal one-qubit 
manipulations. 

VI. TWO-QUBIT GATES 

For the two-qubit gate we consider two ions trapped in the same harmonic well. If the 
frequency difference of the two Raman beams is close to a motional frequency Um and much 
smaller than the qubit level frequency difference, 60 <^ uq, we can neglect terms rotating 
with ujq in Eq. ([1]). Then the effective Hamiltonian of the ions can be approximated by 

H,. = h[n^{r,)\ T 1)(T 1| + l^i(ri)| i l|]e-(^«*+*)e^^'^-^ +h.c. 

+ ;i[fit(r2)| T2)(T2|+fii(r2)U2)a2|]e-^(^°*+<^V^'^-^^+h.c., (11) 

where ri and r2 are the positions of the two ions. We use a shorthand notation for the two- 
qubit operators | t 1| |/2)(/2| = | t l)(t 1| and so on, suppressing identity operators 
/. Two ions in the same harmonic potential well perform normal-mode oscillations around 
their equilibrium position. We choose coordinates such that the axis along which the ion 
conformation aligns coincides with the transport direction (2;-axis). The equilibrium distance 
is determined by the balance of Coulomb repulsion and the restoring forces of the external 
potential and is given hjd= [q^ / {2neom uJqqj^)Y^^, with eo the vacuum permittivity, m the 
mass and ujqom the center-of mass oscillation frequency of the ions j^. We can rewrite the 
ion positions as ri = ro + d/2 e^ + ^ri and r2 = rQ — d/2 ez + 6r2, where Fq coincides with the 
minimum of the harmonic potential, is a unit vector along the 2;-axis and 6rj describes the 
small displacements of ion j around it's equilibrium position. Substituting fq riab + v t, 
we can neglect the overall phase-factor (p + Ak • riab that is common to both ions and, as 
long as it is constant over the time of a gate interaction, does not change the logical phase 
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of the gate. Therefore, in contrast to the one-qubit rotations, no bookkeeping of over the 
trap array is necessary. We include the Doppler-shift Ak ■ \t in 60, define cp = Akzd/2 and 
make the Lamb-Dicke approximation, so that exp[iAk ■ ^r^] ~ 1 + iAk ■ 6rj, which is vahd 
whenever the excursions \Srj\ are small enough that |Ak • 6rj\ ^ 1. We also neglect the 
small differences in Rabi-frequency due to the different positions of the two ions within the 
Gaussian beam (assuming d <^ wq) and set ^si'^j) — ^^s(ro) (■§ e {T? !})• 
Using these definitions and approximations we can rewrite Eq. ffTT]) : 

H,, ~ h[n^ivo)\ T 1)(T 1| +^(ro)| i l|]e-^(''«*+^)(l + zAk.5ri) + C.C. 
+ h [%(ro)| T 2)(T 2| + n^ivo)\ i 2)(i 2|] e-^(^«*-^)(l + zAk ■ 5v,) + c.c. 

(12) 

We can now replace the small displacements 6rj by the corresponding sum of all normal- 
mode displacements, including those in the directions orthogonal to e^. For two ions of 
equal mass the two normal modes along the 2;-axis are an in-phase (center-of-mass, COM) 
oscillation of both ions at frequency ujqom and and out-of-phase (stretch, str) oscillation 
at frequency cUgtr = V^^^com with normal coordinates zqom = {Svi + 5v2) ■ and 

^^str = l/"\/2 {6ri — 6r2) ■ e^. To be specific, we concentrate on detunings 5o ~ tOstr close 
to the stretch mode frequency, so that \6o — Ustrl ^ {^^com, i^str, t^str — ^^com}- In the 
interaction picture, each normal-mode position operator has the same form as the stretch 
mode operator, 

^.tr = j-^(ae-^--* + ate^--*), (13) 

V 2mUstr 

with UstT replaced by the corresponding mode frequency. Here a^{a) is the harmonic oscillator 
mode creation (destruction) operator. Once inserted into Eq. ( JT2i) . all rapidly oscillating 
terms average to zero and are neglected in the following; we keep the near-resonant term, 
which is proportional to e='='('^o-'^str)^ Substituting the normal coordinates, setting rj = Ak • 
Bz^Jh/ (2mc(Jstr) and making the definition 6 = Ugtr — ^0 we get 

^ thrj J2i^ss'e''' at - Al,e-''' a)\ss') {ss'l (14) 

with s,s'e{t,i}, |ss')(ss'| = |sl)(sl| |s'2)(s'2| and 

= —\/2ism{(p) fi-f 
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—V2i sm{ip) 



(15) 



The Hamiltonian H^z generates a time- and internal-state dependent coherent drive 24 1 
resulting from state-dependent dipole forces on the ions. Since the infinitesimal generator of 
the evolution H^z does not commute with itself for different times, and 5 7^ 0, the evolution 
needs to be calculated with either a time-ordered integral approach 25| or another method 
that properly evaluates the acquired state dependent phase. As briefly described in the 
methods section of Ref. ^2^, we can start from the infinitesimal displacement da{t) = 



rjAss'e^^*' generated by H^z and calculate the total coherent displacement 



a{t) = da{t')dt' = r]A,,, f\'^*'dt', 

Jto Jto 



(16) 



and the acquired phase 



Im 
Im 



to 



a*{t')da{t')dt' 



f { t da*{t") dt" I da{t') dt' 

Jtn \Jtn I 



'to \Jto 

to express the evolution Uzz{t) = exp(i$) exp(z[a;(t)a"'' — a*{t)a]) caused by if; 



(17) 



We can 

now restrict our attention to the special case of two Gaussian beams of equal frequency that 
have approximately orthogonal polarization and counter-propagate with ki = — k2 = k (a 
polarization gradient standing wave, see Fig. [2] (b)), with angles 7 and vr + 7, respectively, 
to the transport direction of the ion pair. In the frame of reference of the moving ions, the 
beams are Doppler shifted by Auj = ±|A;| cos(7), respectively, setting the relative detuning 
of the beams to 

5o = 2|A;||t;|cos(7). (18) 

As in the discussion before Eq. (Q we assume that both beams have the same transverse 
mode function given by Eq. and the trajectory is such that the ions traverse the center 
of the beams at t = 0, where they experience the maximum coupling strength fls,m with 
s e {t) i}- The time dependent Raman-coupling strength is then 



= ^^s,mexp[-2(f sin(7)t/t(;o)^ 



(19) 
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We can define the approximate transit duration of tlie ions tlirougli tlie beam r = 
f sin(7)). Since tlie Gaussian envelope is common to both fi^, we can rewrite Eq. 

(USD as 

A,At)=Ks'e-''''\ (20) 
with = As^s'{t = 0)- It is then straightforward to solve the integral for the mode 
displacement in phase space Eq. ( |T6l) 

i5T\ „ I , idr 



t 2 



erf [t/r " ^ j + erf (^t/r + ^ 



(21) 



For the final displacement, after the ions have traversed a distance equal to several beam- 
waists we can evaluate Eq. (!20|) at t = oo and find 



a. 



s'(oo) = e ~ y/iTr]A°^^,r. 



(22) 



This expression shows that the end point of the trajectory can be ass' (oo) 7^ 0. However, 
by choosing 5r/\^ large we can bring the end point exponentially close to zero. This is 
important since we need the motional state of the ions to be disentangled from the internal 
states after the gate is executed. The simplest way to achieve this is to make sure that the 
final displacement for all combinations of internal states is zero, or very close to zero. We 
return to this point and its implications for the gate fidelity below. 

It is also worth noting that from Eq. fl2T]) . ^^^'(oo) is proportional to the Fourier transform 
of the pulse envelope at frequency 6 . This implies that we can make the final displacements 
very small for any choice of envelope that has small Fourier components at 6, and the 
transport gates also work for beam shapes that deviate from Gaussian as long as they are 
reasonably smooth. 

For the logic phases after a complete transit through the beam we need to solve the 
integral 



^ss' = IvAW^m 
= IvAW^m 



oo / rt' 



-co \J—oo 
oo rt' 



dt" e-*"/-'+*^*' dt' 



CO J —CO 



(23) 



A coordinate transformation t" = [u — v)/\/2, t' = {u + v)/\/2 lifts the interdependence of 
the integrals, leading to a closed form solution 



^ss' = IvAW^m 



^"1^" du / e-^'/^'+^^^^ dv 
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= \VAI,? r' " ^ --\V-2J ^ (24) 

where eTfi{z) = —i erf (iz). Figure [3] shows ass'{oo) and both normahzed to their 
respective maximal values as a function of the parameter p = 5t/\^. It gives insight 
into the trade-off necessary to meet the requirement |ass/(oo)| ^ 0. While ass'{oo) decays 
exponentially in p^, decays more slowly with p due to its dependence on the complex 
error function, so that even for very small values of ^^e value of ^ss' remains at 

an appreciable fraction of its maximal value. Note that we could also use two well aligned 
transits through two beam pairs (or the same beam pair). The acquired gate phase would add 
for the two passes, while the displacements can be arranged to be opposite, so the motional 
state returns to the origin. However, for this approach to be viable the timing would have 
to be controlled to a small fraction of 27c/6, which seems imposing, while a sufficiently small 
|ass'(cxD)| turns out to be easily achieved for reasonably large p (see below). The resulting 
phase gate differs by only one-qubit z-rotations from the standard phase gate 

U{^l) = I TT)(TT I + I n)(Ti I + I iT)(iT I + e^*^| |, (25) 

and the total logical phase of the gate can be written as $l = + $^ — ($|j + $||) 
We find 

$^ = -vre-^ erfi (^^^ ^r^'i^l^ - fij V cos(2^). (26) 

This expression shows that the largest modulus of the logical phase is achieved for cos(2y9) = 
±1, corresponding to Ak^d = mr with n integer. Therefore the ion spacing should be a 
multiple of the half-period of the polarization gradient wave along the transport direction z. 
From now on we assume the ion distance is adjusted (by a proper choice of the curvature of 



2 TT erfi ( ^ 



the transported harmonic well) to yield a maximum modulus for the logical phase 26|, |28 |. 
This implies that = = and jA^J^ = \Al^\^ = l/2(fi| - ilif. For a gate equivalent 
to a TT-phase gate we then require $l = — vr, equivalent to 

e-'^^r^^(^^T-fii)Verfi(i|)=l. (27) 

With this condition, and assuming that the trajectory of the ions starts at <y^i{t = — oo) = 0, 
we can derive an expression for the phase space trajectory a-fj(t) that depends only on p 



2^erfi(p) 
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If we assume that the only source of error in the gate is the incomplete return to the 
origin for states that acquire a phase, the fidelity of a gate can be bounded from below by 
the worst-case overlap of the final state with the trap ground state, 

F > mm\J2css' (a,,,(oo)|0)p = min | J] c,,,,e-l"»^'(°°)l'/T = g-l^tiMI' ~ 1 - \au{oo)\^ 

(29) 

where Cggf cirG the amplitudes for state \ss') in an arbitrary two-qubit spin state the gate is 
operating on. Using Eq. (1271) and Eq. (^11) we find 

£= l-F<7r/erfi(p). (30) 

We assume that subsequent gates are implemented after sympathetic re-cooling to the mo- 
tional ground state, preventing coherent addition of errors of this type. 

The trajectories in phase space of the | ||)-state, as described by Eq. (128!) for phase gates 
with e < 10-2 (p = 2.69), e < 10^^ (p = 3.48) and e < {p = 4.11), are shown in Fig. 
m For large \t\, the trajectory spirals near the origin and the number of windings per time 
interval r increases as p is increased. Any desired proximity to the origin can be realized by 
making p large enough. We can therefore also think of p as a gate-adiabaticity parameter, 
because no motional energy is deposited in the system for p — oo. 

In practice it is interesting to note that ^qt/v^ = kwocoti^j) is independent of the 
transport speed and is basically a measure of how many periods of the polarization gradient 
wave are encountered by the ions while they are traveling through the central portion of the 
beam. Unless the beam is very strongly focussed, we have kwo ^ 1, and the ions sample 
many periods of the polarization gradient wave. By choosing 7 close to 90 degrees We can 
reduce the number of periods encountered, but there are drawbacks to such a choice. On 
the one hand, rj oc cos(7) becomes small. As can be seen from Eq. (p7|) a smaller 77 has to be 
compensated somehow, for example, by a higher Rabi rate that requires more laser power, 
could lead to excitation of off- resonant terms in Eq.( fT2l) and cause increased decoherence 
from spontaneous emission |i_6i. In addition, to fulfill the requirement cos(20) = ±1, we must 
restrict the angle 7 to discrete values given by 

Akzd = 2k cos(7„) [g V (2vreomu;coM)] = ^'^^ (31) 

with 7 closest to vr for n = 1. We can still choose different values for 7 by not requiring 
cos(20) = ±1, but at the price of achieving a less than optimal logical phase. 
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To assess the practical feasibility of such a gate, we assume parameters for two ^Be+ ion 
qubits that are achieved in current experiments at NIST, namely an axial trap frequency of 
i^coM = 27r X 4 MHz, a beam waist of Wq = 20 /im, = —1/2Q^, and a wavelength of 313 
nm. We further impose p = 3.48 to ensure that e < 10~^. For n < 10 the necessary detuning 
6 ranges from 29.63 MHz to 768 kHz, too large to reasonably fulfill \6\ <^ ^^com! therefore 
only a more complete theory that takes the faster rotating terms into account in Eq. (fT2|) 
would yield meaningful results. As tabled shows, the transport speed f„ is reasonable for all 
n > 10, and gates with duration r„ on the order of 1 /is can be realized if enough laser power 
is available to produce the Rabi frequencies fij. Gate times are around 5 /is, comparable to 
the present state of the art with phase gates, require 7 to be between 60 and 40 degrees and 
Rabi-frequencies that are close to values that have been realized in previous experiments at 
NIST. 



VII. EFFECTS OF SYMPATHETIC COOLING 



An integral part of QIP with ions in a 



qubit ions with "refrigerator" ions 
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arg e trap array is sympathetic cooling of the 
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32| to remove excess motional energy after 



transport and to reset the motional modes into a well defined initial quantum state. Sym- 
pathetic cooling requires the presence of neighboring ions during gate operations; these were 
not considered in the preceding sections. Although the presence of extra ions complicates 
the description of two qubit gates, it does not change the essential features of the method. 
There are more normal modes of motion in the extended ion configuration and it is more 
complex to properly describe the individual motional amplitudes, especially in a configura- 
tion with ions having considerably different masses. Despite these complications the basic 
gate mechanism is almost unchanged. We consider an ion configuration with a motional 
eigenmode with frequency and normalized eigenvector v. We denote the components 
(amplitudes) of v corresponding to the ion-qubits as Vi and V2 (both are real numbers) 
and substitute these into Eq. (fT2l) . In complete analogy to the calculation following that 
equation, we arrive at the generalized displacement coefficients 

= {vie-''^ + V2e''P) fit 
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All = ivie-''' + V2e'n ^l- (32) 
Calculating the logical phase we obtain the analog to Eq. fl26l) . 

$^ = vre"^ erfi (^^^ 2viV2 r/2(^ - ^ifr'^ cos(2y?). (33) 

Again the logical phase is optimal if the spacing of the qubit ions di2 is such that Akzdi2 = nir 
with n integer. The suitability of a certain mode for a phase gate can be judged by the factor 
V1V2 rf. Both amplitudes should be as large as possible, while rf oc 1/uJv should be not too 
small. At the same time, the refrigerator ions should not have components in v that are too 
small, so that sympathetic cooling is efficient. These requirements are probably best satisfied 
with two ion species of comparable mass in a configuration with refiection symmetry in the 
2;-direction (for example, (g, r, r, q) or (r, g, g, r), where q denotes a qubit and r a refrigerator 
ion). 

The normal modes are then still eigenstates of the reflection operation around z = 0, 
and therefore Vi = ±V2- For "stretch- type" modes with vi = —V2, the optimal coefficients 
fulfill A^^ = A^ = and l^^jP = I^J-fl^ = ^^i(^ — ^i)^; furthermore these modes should 
be better protected from heating because they couple only to field fiuctuations that have an 



appreciable gradient over the ion configuration 331] 



One interesting aspect of sympathetic cooling in the strongly adiabatic regime {p ^ 1) 



was already pointed out by S0rensen and M0lmer [3J]. Due to the small deviation of the 
motional wave packets from the origin and the quick succession of revolutions in phase 
space, the gate becomes more robust to photon recoil and heating of the motion (see also 
|6|]). Therefore such gates can be executed with reasonable fidelity as long as the heating 
rate (in quanta per second) is small compared to the rate of revolution in phase space that 
is set by 6, more precisely, [dn/dt) <^ 1 where is the time to complete one revolution 
in phase space. Gates can even be executed while the ions are re-cooled sympathetically, 
as long as the combined rates of phonon scattering due to heating and cooling are small 
compared to 6. This could be advantageous when very small trap structures are used and 
the ions are in close proximity to electrode surfaces, since the heating rate has been observed 



to scale strongly with the distance to the nearest surface [35|, 136|, 137|, 138| . Sympathetic cooling 
can ensure that the motional state stays in the Lamb-Dicke regime, even if the unchecked 
heating rate would drive the ions out of this regime during a gate operation. 
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VIII. EXTENSIONS 



Hyperfine ground states with an energy splitting that is first-order-insensitive to fluc- 
;uations in the magnetic field have nearly identical AC-Stark shifts, leading to Q-^ ^ 



39 



40 



4l|. As can be seen from Eq. (12^ only a negligible logical phase can be acquired 



in this situation. This problem can be solved in different ways. If single-qubit rotations of 
high fidelity are available, one can transfer out of the field-independent manifold into states 
with 7^ Qi for the duration of the two-qubit gate operation, followed by transferring 
back. The ambient magnetic field fiuctuations typically do not vary appreciably over the 
gate time; therefore phase errors can be sufficiently corrected by spin-echo techniques. The 
extra transfer operations are a disadvantage, but such a scheme has the benefit that qubits 
that are not scheduled for gate operations are not as readily dephased by Stark shifts caused 
by residual stray light from two-qubit gate beams, which are likely to be the most powerful 
laser beams used in a working QIP device. 

Another approach is to use the gate first described by S0rensen and M0lmer 421] and 
Solano et al. [43l| . The gate mechanism is basically the same as described above, but in a 
rotated basis |41[ ^ and can be adapted as a transport gate in the more robust way described 
in Fig. 1(b) of j4l|. Here, two counter-propagating beams are offset from the primary 
beam by ±ujq as shown in Fig. [51 For the prop er transport velocity, this would set up a 
situation equivalent to the one described in 4l|. The frequency differences resulting from 
combining these three Raman beams, including the appropriate Doppler-shifts due to the 
transport, are at Uq ± Ugtr + creating slow rotating terms near the blue uq + Ustr and 
red ujq — cUstr sideband transitions in the Hamiltonian ([T]). This scheme has the advantage 
that it works for qubits with field-insensitive qubit transition frequencies. However, this 
comes with some disadvantages. One is the added technical difficulty inherent in creating 
the frequencies of the required beams. For example, the scheme depicted in Fig. [5] would 
require acousto-optical modulation at u;o/2. The smallest value of ujq for typically considered 
hyperfine-state ion-qubits is 2tt x 1.25 GHz in ^Be"*". Even in this case we would need 
two specialized modulators at approximately 600 MHz, adding technical complexity. For 
hyperfine-splittings exceeding 2 GHz the modulation would currently require several acousto- 
optic and/or electro-optic modulators, or a means to offset-phase-lock the respective laser 
beams. Another disadvantage is the relative proximity of the qubit carrier transition at the 
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difference frequency ujq. To keep perturbations from carrier excitation small, we have to use 
gate durations much longer than 271/^7811- For the Z-phase gate described in the previous 
sections, the carrier is essentially detuned by ujq, so gate durations comparable to or shorter 
than 27r/co'coM could be possible. 

Another possible extension consists of fine tuning the phase-space trajectories with trans- 
port patterns featuring nonuniform velocities. Such schemes could lead to better suppression 
of off-resonant terms and possibly to higher gate speeds. This would probably require nu- 
merical modeling, while the simple uniform-velocity transport through Gaussian beams has 
the appeal that all relevant quantities can be expressed analytically. 

More generally the transport gate mechanism can be implemented with any standing wave 
field that imparts an internal state-dependent force. For example, close proximity of the 
ions to the trap electrode surfaces should make it possible to set up strong magnetic field 
gradients to induce state-dependent forces. We could construct a magnetic "washboard" 
resulting from parallel current-carrying wires with alternating current direction, or strips 
of permanent magnetic material with alternating polarization oriented perpendicular to the 
transport direction. As long as the qubit information is contained in states with different 
energy dependence (slope) as a function of the external field, the field gradients of the 
washboard produces state-dependent forces. By controlling the transport velocity of the 
ions we could make them experience these forces in near resonance with a motional mode 
frequency. Once this situation is established, the gate mechanism is analogous to what has 
been described in the previous sections, where the strength and envelope of the alternating 
force would be given by the particular values of the currents used or the strength of the 
magnetized strips. Thus a range of envelope functions could be achieved. For example, 
assume a washboard made of FePt stripes (or domains written into a continuous film) with 
an alternating magnetization on the order of 800 kA/m, and a period length of dm — 20 fim 



could be produced (similar to those of Ref. [4J]). As a ^Be"*" ion is transported at a height 



of 17.5 /xm over such a structure at vw = 80 /xm//is, this would create a time- dependent 
magnetic field with an amplitude of 



Btot{t) = ^{Bo - B^cos{uj^t)y + {B^sm{uj^t)y - Bo{l + 1/A{B^/BqY) - 5^cos(cu^t), 

(34) 



with R„ ~ 20 G, 5n ~ 120 G 



45| . This choice of Bq ensures first-order field independence to 
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a qubit encoded in the \F = l,mF = 0) and |F = 2, m^? = 1) states in all regions of the array 



that are sufficiently far from the permanent magnetic elements j40|. To lowest order we have 
an oscillating field at = 2TiVw/dm = 4 MHz that is superimposed on a static average field 
that is larger than Bq by 5^/ (4i?o) — 0.837 G. We assume that good single-qubit rotations 
are available to transfer the qubit information into the \F = 2, nip = —2)/\F = 2, nip = 2) 
"stretched" states of the hyperfine manifold. Then the undulating field of the washboard 
would lead to a gate Rabi frequency of 

f]„ = ^^^/^ ^ 27r X 73.5 kHz, (35) 



with the Bohr magneton and Zq = yh/{AniuJcoM) the extent of the ground-state wave- 
function of the COM mode of two beryllium ions. This Rabi frequency would correspond 
to a gate duration Tm = ir/^m — 6.8 fis. At the same time, the slightly higher average field 
leads to effective single-qubit rotations of about 27rx 15.9 that commute with the phase gate. 
These unwanted single-qubit phases could possibly be eliminated by a more sophisticated 
(symmetrical) geometry of the elements producing the oscillating field, or by dividing the 



gate into two interactions with was^ 



enclosed in a spin-echo sequence [26 1. 



iboards that each yield a logical phase of 7r/4 and are 



Gate speeds are limited by the condition that the spins should follow the magnetic field 
adiabatically. For this to happen, the normalized rate of change of the field \{dB/dt)\/ Bq 
has to be much smaller than the smallest Larmor frequency along the trajectory ujmin = 
min(7<|B|), with 7< the smallest gyromagnetic ratio of the states in question and B the 
position-dependent magnetic field. The speed of variation of B is governed by the trap 
frequency, therefore, at least for a gate utilizing adiabatic following, the gate duration is 
limited by the trap frequency. The "Lamb-Dicke" parameter that is relevant for this gate is 
rjm = {2n/djn) zq. In contrast to gates mediated by laser fields with wavelengths prescribed 
by the internal states of the ions, this parameter can be scaled more freely by changing 
the period length dm of the microfabricated pattern. Of course the ions have to be moved 
closer to the surface to create a sufficient modulation of the field along their trajectory. This 
might lead to an increase in anomalous heating, but at the same time the closer proximity 
of the ions allows use of smaller trap electrodes and scaling the dimensions of the trap array. 
For gates mediated by light, the need to illuminate ions by laser beams with waists that 
are most likely limited to sizes of several micrometers from geometrical constraints dictates 
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a minimum distance of the ions to the trap electrodes. This distance in turn governs the 
minimal electrode dimensions in a trap array. If light is not necessary to mediate the gates 
we could use trap structures beyond this limit, possibly with stronger confinement of the 
ions and faster gates. 

An additional, but important advantage of such gates is the absence of spontaneous 
emission that is currently thought to pose the most severe limitation on achievable fidelities 



for Raman transition gates implemented with laser beams 



6|- 



IX. SUMMARY AND CONCLUSIONS 



Miniaturization of ion traps and scaling up ion qubit numbers to beyond 10 is now under 
way in several laboratories. Recent efforts have concentrated on developing technologies for 
trap arrays compatible with this goal. One of the next steps is to simplify the optics necessary 
for scale-up as much as possible. In this paper we have outlined a few possible techniques 
towards this goal. Controlled transport of ions is utilized for parallel implementation of 
quantum logic gates with relaxed requirements for temporal and spatial control of the laser 
beams. After sketching the basic features of such an architecture based on a multi-zone trap 
array, we analyzed one-qubit rotations and a universal two-qubit gate and showed that such 
gates can be implemented with existing technology. 

We also discussed possible extensions of the two-qubit gate mechanism. In particular we 
sketched one approach that utilizes periodic magnetic field patterns in combination with 
ion transport to exert state-dependent forces on the ions. This approach might enable more 
flexible scaling of future trap arrays and could remove some of the limitations to fidelity and 
gate speed posed by two-qubit gates mediated by laser beams. 
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TABLE I: Gate parameters for an axial trap frequency of cJcoM = 27rx 4 MHz, beam waist of wq = 
20 fim, fi-f = —1/2 ^l^, wavelength of 313 nm (^Be+) and p = 3.48 {e < 10~^). All parameters, 
including Rabi-frequencies Q^, are within experimental reach. 
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FIG. 1: (Color online) Basic steps for the proposed architecture: (a) The ions carrying the quantum 
information are arranged into the desired spatial configuration in the trap array while the laser 
beams are switched off. (b) All laser-beam assisted operations scheduled after the prearrangement 
in (a) are carried out. This includes one-qubit rotations, two-qubit gates and measurements. 
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FIG. 2: (Color online) Beam parameters for one- and two-qubit gates: (a) Orientation of the 
wavevectors ki and k2 relative to the transport direction for co-propagating beams as used in 
the one-qubit rotations, (b) Orientation of the wavevectors ki and k2 relative to the transport 
direction for counter-propagating beams as used in the two-qubit phase gate, (c) The starting 
point of the ion trajectory is a distance D from the center of the Gaussian beam profile at z = 
with waist size wq. 
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FIG. 3: (Color online) Normalized modulus of the final coherent state amplitude a„ 
|qss'(oo)|/|qss/(0)| (dashed) and normalized phase = ^ss' / ^^^[^ss'] (solid) as a function 
P- 
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FIG. 4: Phase space trajectories of the state | Ti) for (a) e < lO'^ {p = 2.69), (b) e < 10~^ 
{p = 3.48) and (c) e < 10~^ {p = 4.11). The arrows indicate the trajectory directions. The 
trajectories are fully determined (up to their orientation in phase space) by p and the requirement 
that the logical phase is $l = — tt. 
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FIG. 5: Possible configurator! for the transport version of a S0rensen and M0lmer gate. The two 
double-pass acousto-optic modulators (AOMs) shift the retro-reflected beams by ±a;o, respectively. 
Together with the Doppler shift due to the transport, the ion experiences beat notes close to its 
blue and red sidebands on the stretch mode that implement the gate. 
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